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Abstract 

We perform fully coupled two-dimensional numerical simulations of plane channel helium II counter¬ 
flows with vortex-line density typical of experiments. The main features of our approach are the inclusion 
of the back reaction of the superfluid vortices on the normal fluid and the presence of solid boundaries. 
Despite the reduced dimensionality, our model is realistic enough to reproduce vortex density distributions 
across the channel recently calculated in three-dimensions. We focus on the coarse-grained superfluid and 
normal fluid velocity profiles, recovering the normal fluid profile recently observed employing a technique 
based on laser-induced fluorescence of metastable helium molecules. 
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I. INTRODUCTION 


Three-dimensional homogeneous isotropic turbulence is the benchmark of turbulence research. 
Recent papers^! have compared the properties of homogeneous isotropic turbulence in ordinary 
(classical) fluids and in liquid helium near absolute zero, and found remarkable similarities. In 
particular, experiments have revealed that the temporal decay of vorticily- is the same, and that the 
energy spectrum (which represents the distribution of kinetic energy over the length scales) obeys 
the same classical Kolmogorov scaling at sufficiently large length scales,-- in agreement with 
theoretical^ and numerical studies.-^— These results are surprising, because the low temperature 
phase of 4 He (hereafter referred to simply as helium II), is quite different from an ordinary fluid.— 
It is well-known, in fact, that helium II consists of two interpenetrating fluid components: a 
viscous normal fluid (whose vorticity is unconstrained) and an inviscid superfluid (whose vorticity 
is confined to vortex line singularities of fixed circulation h/m where h is Planck’s constant and 
m the mass of one helium atom). 

Despite its importance, isotropic homogeneous turbulence is an idealization which neglects the 
role of boundaries (for example, vorticity is generated at the walls of a channel). In this report we 
are concerned with superfluid turbulence along channels or pipes. Such flows are neither homoge¬ 
neous (because the boundary conditions are likely to induce non-uniform profiles) nor isotropic 
(because of the direction of the flow). The prototype channel problem of the helium literature is 
thermal counterflow.—The typical experimental set-up consists of a channel which is closed at 
one end, and is open to the helium bath at the other end. At the closed end, a resistor dissipates a 
known heat flux which is carried away by the normal fluid; to conserve mass, superfluid flows in 
the opposite direction towards the resistor; the resulting velocity difference between the two fluids 
is proportional to the applied heat flux. If this heat flux is larger than a small critical value, the 
superfluid component becomes turbulent, forming a disordered tangle of quantised vortex lines 
(superfluid turbulence). The intensity of the vortex tangle is usually characterized by its vortex 
line density L (length of quantized vortex lines per unit volume), which can be determined by 
measuring the attenuation of second sound as a function of the applied heat flux. 

The questions which we address in this work is simple but fundamental: what are the profiles 
of the normal fluid, of the superfluid, and of the vortex density across the counterflow channel? 

This question motivates current experimental attempts to directly visualize the flow of he¬ 
lium II. Two new visualization methods stand out. Particle Tracking Velocimetry (PTV) of hydro- 
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gen and/or deuterium flakes^r— has been used to image individual quantum vortex reconnections^ 
and to determine the velocity and acceleration statistics of the turbulent supcrfluid.-^^ Laser- 
induced fluorescence of metastable helium moleculehas directly imaged the profile of the 
normal component, addressing the issue of whether, at sufficiently large heat currents, the normal 
fluid flow undergoes a laminar-turbulent transition.— 

Until now, the question of the profiles of normal fluid, superfluid and vortex line density has 
been unanswered. On first thoughts, in analogy with a classical viscous fluid (which obeys the 
Navier-Stokes equation with no-slip boundary conditions), the normal fluid component should 
have a parabolic Poiseuille profile across the channel; similarly, in analogy with a classical inviscid 
fluid (which obeys the Euler equation and, unimpeded by viscosity, can slip along the channel’s 
walls), the superfluid component should have a uniform profile, hence the vortex line density 
should be non polarized and, eventually, uniform. On second thoughts, the said profiles cannot 
be correct: the varying mismatch between the superfluid and normal fluid velocities across the 
channel would induce a large non-uniform mutual friction^ which would modify these profiles. 
To appreciate the mathematical difficulty of the problem, notice that not only are the two fluid 
components coupled, (the normal fluid affects the superfluid and viceversa), but the coupling term 
between the two fluids is nonlinear: the mutual friction is proportional to the velocity difference 
between normal fluid and superfluid, times the vortex line density, which is a nonlinear function 
of this velocity difference. 

Unfortunately, most numerical simulations of superfluid turbulence in the literature have de¬ 
termined the superfluid vortex tangle in the presence of a prescribed normal fluid, without tak¬ 
ing into account the back reaction of the vortex lines on the normal fluid. Various models of 
the imposed normal fluid have been studied: uniform,— parabolic,— — Hagen-Poiseuille and 
tail-flattened flows,— vortex tubes,— ABC flows,— frozen normal fluid vortex tangles^, random 
waves^, time-frozen snapshots of the turbulent solution of the Navier-Stokes equation s 35 ' 37 ’ 39 and 
time-dependent homogeneous and isotropic turbulent solutions of linearly forced Navier-Stokes 
equations.— Moreover, most calculations were performed in open or periodic domains, avoiding 
the difficulty of the boundary. Other works have determined the effects of a prescribed superfluid 
tangle on the normal fluid,— failing again to fully model the coupling of superfluid and normal 
fluid. Because of the computational complexity and cost involved, fully coupled calculations have 
been attempted only for simple configurations, such as single, isolated vortex lines^ or rings,— or 
for decaying tangles in open geometry^ and periodic domains 4^ 
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The model which we present here is fully coupled (the normal fluid affects the superfluid and 
viceversa via a nonlinear mutual friction term) and includes boundaries. To cope with the com¬ 
putational difficulty, our model is two-dimensional rather than three-dimensional: vortex loops 
in a three-dimensional channel are thus replaced by vortex points in a two-dimensional channel. 
Despite the simplified dimensionality, our model captures the nonlinearity of the problem, which, 
we think, is the key ingredient to determine flow profiles in actual channels. 

The outline of the paper is the following. In Section [II] we describe the two-dimensional model 
which we use and the details of the numerical algorithm. Section [III] focuses on the results and in 
Section [TV] we critically discuss to what extent our two-dimensional model is capable of grasp¬ 
ing the most relevant vortex dynamics occurring in helium II counterflows. Finally, Section [VI 
summarizes the conclusions. 

II. MODEL 

A. The counterflow channel 

We consider an infinite two-dimensional channel of width D. Let x and y be respectively the 
directions along and across the channel with walls at y = ±D/2 and periodic boundary conditions 
imposed at x = 0 and x = L x . The average normal fluid and superfluid flows are respectively in 
the negative and positive x direction. 

The superfluid vortices are modelled as N vortex-points of circulation T ;; and position r j(t) = 
(xj(t), yj(t)), where j = 1, - ■■ N and t is time. 

Half the vortices have positive circulation Yq = k and half have negative circulation T ; = — k, 
where n = 10 _3 cm 2 /s is the quantum of circulation in superfluid 4 He. 

To make connection with experiments we interpret n = N/ ( DL X ) (average number of vortex- 
points per unit area) as the two-dimensional analogue of the three-dimensional vortex-line density 
L, and relate L to the channel-averaged normal fluid longitudinal velocity (u n ) via the relation-^ 


L 1/2 D = 1.03^0—(u n )h D - 1.48/3, 


( 1 ) 


Ps 

where (u n ) is related to the applied heat flux q via 



( 2 ) 
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where T is the absolute temperature, S the specific entropy, and p n , p s and p the normal fluid, 
superfluid and total helium II densities, respectively, where p = p n + p s . The coefficients 70 
and /? in Eq. (Q]) have been determined experimentally by Tough and collaborators,— — while ho 
represents the channel’s hydraulic diameter. 

In the absence of vortices, the counterflow condition of zero net mass flow 

Pn(u n ) + PsV £ s Xt = 0, (3) 

determines the uniform superflow v e s xt in the opposite direction with respect to the normal fluid. 
Notice that Eq. ([[]) coincides with Vinen’s equation^ describing the evolution of the vortex-line 
density L modified in order to take into account the presence of solid boundaries and that the 
average intervortex distance i is defined by the relation = L -1 / 2 . 

B. The superfluid vortices 

The vortex points move according to22 


-jj- = v s (r j, t) + a s' x (v n (r i , t) - v a (r t)) 

+a' (v n (r,-, t) — v s (rj,t)) (4) 

where s' is the unit vector along vortex j (in the positive or negative z direction depending 
on whether Tj is positive or negative), a and a' are temperature dependent mutual friction 
coefficients^, v n (r j,t) is the normal fluid velocity at position r, ; the superfluid velocity at position 
rj is decomposed as 

v s (r j,t) = v e s xt (t) + v si (r j: t), (5) 

where v( ,T ' ( t) is the uniform (potential) superfluid flow which enforces the counterflow condition 
of no net mass flow and v si (r j, t ) is the superfluid velocity field induced by all the N vortex-points 

at r r 


v 8i (r j ,t)= ^ v s?;,fc( r j! t). ( 6 ) 

k=l...N 

The integration in time of Eq. © is performed employing the second-order Adams-Bashfort 
temporal advancement scheme. 
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To determine the superfluid velocity field induced by the k -th vortex v.^/Jx. t) we employ 
a complex-potential-based formulation enforcing the boundary condition that, at each wall, the 
superfluid has zero velocity component in the wall-normal direction. 

The complex potential can be derived using conformal mapping^ or, equivalently, using (for 
each vortex) an infinite number of images with respect to the channel walls,— leading to the fol¬ 
lowing expression 


F k (z,t ) 




h 


2 7T777, 


log 


sinh [^{z - z k {t))\ 
sinh [JL.( z -zi(t))] 


(7) 


where z k (t) = x k [t) + iy k (t) is the complex number associated to r k (t). The corresponding 
superfluid velocity v sijk (z, t ) = (y x si k , v y si k ) is obtained from the complex potential in the usual 
way as 


J si,k 


• v 

- IV * u = 


dF k (z,t ) 


( 8 ) 


si,k dz 

The uniform superfluid velocity v“ 4 (f) = (■ Vg Xt (t ), 0) in Eq. © is instead obtained by enforc¬ 
ing at each timestep the counterflow condition of no net mass flow taking into account the presence 
of vortices, i.e 


Pn{u n ) + Ps ({Usi){t) + uf*(t)) = 0 . (9) 

where v si = (u sl , v si ) to ease notation. 

To model the creation and the destruction of vortices (mechanisms intrinsically three-dimensional) 
within our two-dimensional model, we proceed as follows. When the distance between two vortex 
points of opposite circulation becomes smaller than a critical value e\, we perform a ’’numerical 
vortex reconnection” and remove these vortex points; similarly, when the distance between a 
vortex point and a boundary is less than 62 = 0.5ei, we remove this vortex point (the vortex of 
opposite circulation being the nearest image vortex beyond the wall). To maintain a steady state, 
when a vortex point is removed, a new vortex point of the same circulation is re-inserted into 
the channel in a random position. In order to assess the dependence of the numerical results on 
the value of e 1; we have performed numerical simulations varying the value of e\ by two orders 
of magnitude: we find that the results are identical. This reconnection model, corresponding 
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three-dimensionally to the vortex filament method of Schwarz,— correctly describes the fate of 
two very near antiparallel vortices (as confirmed by past Gross-Pitaevskii numerical studies^) and 
avoids the generation of infinitesimal length scales which would trigger numerical instabilities. 
In order to estimate the impact of this re-insertion procedure on the numerical results, another 
two-dimensional renucleation model has also been explored in the present study: the vortices are 
re-inserted with the same wall-normal coordinate y with which they have been removed and a 
random streamwise x coordinate. The results obtained are quasi-identical to the ones obtained 
with the random re-insertion model, concluding that the numerical results presented in Section 
IIIIB1 referring to the random re-insertion model, are not an artificial outcome of the reconnection 
procedure. 

C. The normal fluid 

Typical experimental values of pressure and temperature variations along counterflow channels 
allow us to assume that both superfluid and normal fluid flows are incompressible and isoentropic, 
i.e. p . [),,. p s . S arc constant. Furthermore, assuming negligible the variations of the normal fluid 
dynamic viscosity r] n and of the thermal conductivity A across the channel and neglecting quadratic 
or higher-order terms in spatial gradients of velocity and thermodynamics variables, the resulting 
incompressible and isoentropic equations of motion of the normal fluid are the following:^^ 


<9V„, 

dt 


+ (v n • V) v n 


—-Vp — —SWT + r' n V 2 v„ 
P Pn 

_ 7 tV (v„ - v s ) 2 + —F ns 
2 P pn 


( 10 ) 


V ■ v n = 0 (11) 

where v n = r\ n fp n is the normal fluid kinematic viscosity, and the mutual friction force F n , s is 
determined by the averaging procedure described in Section HlDl 

The normal fluid velocity field v„ is decomposed in two solenoidal fields: 


v 


n 


+ V . 

n ' n 


( 12 ) 
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The first field = (u y , v%) = (— V n0 [1 — (2 y/D ) 2 }, 0) is the Poiseuille flow which would exist in 
absence of superfluid vorticity at constant heat flux q supplied by the heater. The second velocity 
field v' ( = (?//.,, v' n ) accounts for the back reaction of the superfluid vortex-lines on the normal 
fluid. To calculate v' ; we employ the vorticity-stream function formulation, according to which 
we define the stream function T' and vorticity field u' n as follows: 


v 


r 

n 


chk'\ 
dy dx ) 


(13) 


= (V X V^) • Z , 


(14) 


where z is the unit vector in the z direction. The definition of T', Eq. (fl3l) . directly ensures that 
v' n is solenoidal, Eq. (fill) , while the Navier-Stokes equations (fTOl) are equivalent to the following 
two scalar equations: 

W = -J n (15) 


dt 


+ 



du'n 

dx 


(K>' ( duj' n 
dx \ dy 



V72 / | 

V n \7 LO n -\ - 

Pn 


dpy 

dx 



(16) 


where F ns = ^ F x , F y ^J. 

The evolution equation (fT6l) for the normal vorticity c J n is discretized in space employing 
second-order finite differences and its temporal integration is accomplished using the second- 
order Adams-Bashfort numerical scheme. The Poisson equation (fl5l) is instead solved in a mixed 
(k x , y) space, employing a Fourier-spectral discretization in the periodic /-direction and second- 
order finite differences in the wall-normal direction y. The boundary conditions on T' and u' n are 
deduced by imposing no-slip boundary conditions on the viscous normal fluid velocity field. 


D. The mutual friction 

The mutual friction force F ns accounts for the momentum exchange between the normal fluid 
and the superfluid in presence of the quantized vortex-lines which act as scattering centres for 
the elementary excitations constituting the normal component.— This exchange takes place at very 
small length-scales, less than the average intervortex distance l, beyond the practical numerical 


8 












resolution and, at some temperatures, the hydrodynamical description of the normal fluid. To make 
progress, we employ the coarse-grained theoretical framework elaborated by Hall and Vinen^ 
according to which, at lengthscales larger than t, the mutual friction forcing assumes the following 
expression 

F ns = ap s uj s x [u3 s x (v n - v s )] + ol p s u: s x (v n - v s ) , (17) 

where ~ symbols indicate coarse-grained averaged quantities. 

We distinguish between the (Ax, Ay) grid on which the normal fluid velocity v n is numeri¬ 
cally determined, and the coarser (AX, AY) grid on which we define the mutual friction F ns . In 
principle, we would like to have AX and AY /$> l corresponding to the Hall-Vinen limit; in 
practice, we use AX and AY > l due to computational constraints. To prevent rapid fluctuations 
of the friction at small length-scales, we smooth the vortex distribution using the Gaussian kernel 
©j(r) associated to each vortex j according to the following expression 


r — r. 


0 j( r ) = 2£2 


(18) 


Lx D/2 


where V) = 



i r ~ r j i 2 

2< 2 dxdy. Hence, on the basis of Eq. (flTl) . the mutual friction force F 


i p,q 


0 -D/2 


averaged on the coarse grid-cell ( p , q) is given by the following expression 

Fns = -a P s (v*« - 


+a'p s fi M zx(^-v«) 


(19) 


where 


£M = E 


j=l...N 


AXAY 


Qj(r)dr 


(p.?) 


( 20 ) 


n p ’ q = 


j=l:.N 


AXAY 


Qj(r)dr 


(P.9) 


( 21 ) 


Tj = ±k and the symbol jj denotes the integral over the coarse grid-cell (p, q). Physically, 

(p,q) 

L p,q corresponds to the coarse-grained vortex-line density while Q p,q coincides with the coarse¬ 
grained superfluid vorticity. Finally, we average F p / q over the short time interval T ns = AX/v e s xt . 
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the average time interval during which a vortex-point moves from a coarse grid-cell to the neigh¬ 
bouring (cfr. Eq. ©). 

The interpolation of F ns on the finer grid (Ax, Ay) is performed via a two-dimensional bi¬ 
cubic convolution kernel^ whose order of accuracy is between linear interpolation and cubic 
splines orders of accuracy. The structure of the fine and coarse grids on a particular portion of 
the computational domain is illustrated in Fig. Q] while in Fig. [2] we report a two-dimensional 
color plot of the longitudinal component of the mutual friction force F x interpolated on the fine 
grid, on the same domain as Fig. Q] the smoothing effect of the Gaussian kernel combined with the 
interpolating scheme emerges clearly, if compared to the ideally 5-shaped nature of F ns centered 
on the vortex-points displayed in Fig. Q] Furthermore, it is worth emphasizing that the employ¬ 
ment of Eq. (fT9l) for the computation of the mutual friction force F ns , ensures a smooth transition 
when the vortex-points cross coarse grid-cell boundaries. 


III. NUMERICAL SIMULATIONS 


A. Parameters 


We chose the parameters of the numerical simulations in order to be able to make at least 
qualitative comparisons with experiments. As a reference, we select the experimental counterflow 
studies performed by Tough and collaborators on both high aspect-ratio rectangular cross-section 
channels,— which represent the closest real experimental settings to our idealized plane channel, 
and cylindrical capillary tube a —More in detail, we set the width of the channel D = 9.1 x 
10“ 3 cm, corresponding to tube R4 in Ref. [190, and = 25. The consequent Reynolds 

number of the normal fluid flow calculated via Eq. ([]]) is Re = 206, far below the critical Reynolds 
number for the onset of classical turbulent channel flows Re c ~ 5772.— As a consequence, on the 
basis also of past experimental i nvesti gation s ,—— we reckon that in our numerical experiment the 
flow of the normal fluid is still laminar. 

The complete list of parameters employed in our simulation and the subsequent physical rele¬ 
vant quantities are reported in Table [D expressed in terms of the following units of length, velocity 
and time, respectively: 5 C = D/2 = 4.55 x 10 _3 cm, u c = n/{2i t 5 c ) = 3.49 x 10 _2 cm/s, 
t c = 8 c /u c = 0.13 s. Hereafter all the quantities which we mention are dimensionless, unless oth¬ 
erwise stated. The constant V n0 determining u/ is computed imposing, without any loss of gener- 
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Figure 1: (Color online). The structure of the fine (blue solid lines) and coarse grids (green solid 
lines) are illustrated on a particular portion of the computational domain, together with positive 
and negative vortices indicated with empty red and filled black circles, respectively. 
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D 

2 

V n0 

553.6 

Lx 

6 

T 

1.7 K 

N 

1876 

Ps/ Pn 

3.373 

n 

156.3 

ei 

2.5 x 10" 3 

£ 

0.08 

A t v 

7.5 x 10~ 6 


Table I: Numerical parameters employed in the simulations 
and subsequent physical relevant quantities in dimensionless units 
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Figure 2: (Color online). Two-dimensional color plot of the longitudinal component of the 
mutual friction force F x /p n (see Eq. (fl6l) '). interpolated on the fine grid on the same domain as 
Fig. m The axes of the plot are rescaled employing the scaling units defined in Section [Til Al 


ality, that the whole normal fluid flow rate is supplied by the Poiseuille field v£, i.e. (u n ) = (?./('), 
implying (u' n ) = 0. In the spirit of the coarse-grained description illustrated in Section HlDl we 
define a coarse and a fine grid characterized by numbers of grid-points and spacings listed in Table 
HI! satisfying the condition AX, AY > £ > Ax, Ay. 

The coupled calculation of vortex motions and v n entails the simultaneous existence of two dif¬ 
ferent timestep stability criteria, one for each motion. Concerning the evolution equation (fl6l) for 
u' n , the constraint is set by the normal fluid viscosity,— leading to the restriction At n < (Ax) 2 /u. 
Regarding the motion of the superfluid vortices, consistently with the numerical reconnection pro- 
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Fine grid 

Coarse grid 

n x 

192 

N x 

48 

n y 

64 

Ny 

16 

Ax 

3.125 x 10 -2 

AX 

0.125 

Ay 

3.125 x 1CT 2 

AY 

0.125 


Table II: Number of grid-points and spacings in dimensionless units 
of the grids employed in the numerical simulations 


cedure illustrated in Section HlBl the integration timestep A t v for Eq. © must satisfy the condition 
At,; < f ] /(2V ei ), where V ei is the velocity of a pair of anti-vortices along their separation vector 
when separated by a distance equal to e\. This constraint on At,, prevents from the generation 
of unphysical small-scale periodic motions ( e.g . vortex-pairs multiple crossings). The value of 
At„ employed in our simulation is reported in Table U and the viscous constraint allows us to set 
A t n = 2 A t v , implying that vortex motions alternately take place with frozen normal fluid. 


B. Results 

1. Steady-state regime 

The aim of our numerical simulations is to determine the spatial distributions of positive and 
negative vortices and the normal fluid and superfluid velocity profiles across the channel in the 
steady-state regime. To stress that these distributions and profiles are meant to be coarse-grained 
over channel stripes of size AY, we use the “ symbols. Fig. [3] illustrates the initial conditions 
of a typical simulation. Fig. [3] (top) shows the initial random spatial distribution of the vortices, 
corresponding to the coarse-grained vortex density profiles n(y) shown in Fig. [3] (middle). In 
Fig. [3] (bottom) the initial parabolic Poiseuille profile for u n and the flat profile for u s are re¬ 
ported. After a transient interval whose characteristics will be addressed in section IIIIB 21 the 
system reaches the statistically-steady-state described in Fig. [4] As expected, the steady-state 
regime is achieved after a time interval Tj rs l) 1 jv. The most important feature is the shape of 
the coarse-grained profile of the normal fluid velocity u n reported in Fig. @] (bottom), which is 
slightly flattened in the near-wall region and sharpened in the central region with respect to the 
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Poiseuille profile. These characteristics have recently been observed experimentally by means of 
laser-induced fluorescence^ in the same counterflow regime (turbulent superfluid, laminar normal 
fluid). In the experiment, the flattening of the profile is more pronounced, but we reckon that this 
difference is due, at least partially, to a larger superfluid turbulent intensity in the experimental 
setting (40 < L l / 2 D < 70 against L X ^ 2 D ~ 25 in our simulation). 

The other key feature which emerges from the numerical simulation is the polarization of the 
superfluid vortex distribution, which can be qualitatively observed in the snapshot of the steady- 
state vortex configuration, Fig. 0 ] (top). To investigate quantitatively this aspect, we introduce the 
coarse-grained polarization vector p(y) defined by“ 


v a (y) n + (y)-n (y) „ 

- - - ^ 

ttn(y) n+(y) +n~(y) 


( 22 ) 


Note that p (y) = 0 when quantum turbulence is uniformly distributed all over the channel (as, 
for instance, at t = 0 in our numerical simulations, see Fig. [3] (top) and (middle) ). The steady- 
state profile of the polarization magnitude p(y) is reported in Fig. 0] (middle) together with the 
positive and negative vortex density profiles, n + (y) and n~(y) respectively. This polarized pattern 
directly arises from the vortex-points equations of motion ©, where the friction term containing 
a depends on the polarity of vortex. 

This polarization of the vortex configuration, which, we stress, is not complete, i.e. \p(y)\ < 1, 
generates a parabolic coarse-grained superfluid velocity profile u s (y) ~ y 2 which is reported in 
Fig. ED (bottom). This process, i.e. the superfluid polarization induced by a normal fluid shear 
generating a superfluid velocity pattern which mimics the normal fluid one, confirms past analyti¬ 
cal results obtained via simple models^ and backs numerically observed normal fluid-superfluid 
velocity matching and vorticity locking 

It is interesting to notice that our model, although being two-dimensional, recovers the total 
vortex density profile n(y) computed very recently via three-dimensional numerical simulations 
of helium II channel counterflows with prescribed Poiseuille normal flow.—— On the contrary, the 
vortex density profile n(y) calculated in this work is significantly different from the ones computed 
in past two-dimensional simulations with prescribed Poiseuille normal flow, where the density is 
approximately uniform across the channel.-^^ 
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2. Transient inten’al 


The main results of our investigations have been outlined in the previous section. Before we 
finish, it is instructive to describe how the vortices and the normal fluid adjust to each other reach¬ 
ing a steady-state, starting from our arbitrary initial condition: this exercise helps to understand 
the physics of the coupling of vortices and normal fluid. 

The evolution to the steady-state can be understood using the coarse-grained profile of the 
longitudinal component of the mutual friction force F x , reported in Fig. [ 6 ] The expression of F ‘ l 
at a first order of accuracy according to Eq. (fl9l) . is 

F X (y) ^ —ap s Kn(y) ( u n (y ) - u s (y)) (23) 

At t = 0, F‘ is stronger in the central region of the channel, flattening the profile of the normal 
fluid at time t x ~ 6.8 x 10~ 3 7j very close to the initial configuration, as illustrated in Fig. [5] 
(bottom). At times t ~ ti, the superfluid polarization is only partial, see Fig. [5] (top), generating 
a less pronounced superfluid velocity profile u s (y) (Fig. [5] (bottom)). The resulting longitudinal 
component of the mutual friction force at t ~ ti is therefore more uniform across the channel with 
respect to t = 0, as illustrated in Fig. [ 6 J This allows the normal fluid to regain a quasi parabolic 
profile in the subsequent time interval (u n is approximately parabolic at t ~ 0.25 Tf). Finally, at 
t = Tf, the flow reaches a self-consistent dynamical equilibrium determined by (a) the vortex- 
density and velocity profiles reported in Fig. [4] (middle) and (bottom) and (b) the longitudinal 
component of the mutual friction force illustrated in Fig. ( 6 [ characterized by peak values in the 
near-wall region. 


IV. DISCUSSION 

The aim of the present section is to (a) describe the idealized three-dimensional dynamics 
which we reckon corresponds to the two-dimensional vortex-points motion illustrated in Section 
fill Bl and (b) critically discuss to what extent this idealized three-dimensional motion is capable of 
grasping the most relevant vortex-tangle dynamics occurring in helium IIT-I counterflows. These 
two issues will be addressed in Sections [IV Al and llYBl respectively. 
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Figure 3: (Color online), (top): vortex distribution at t = 0, red empty (black filled) circels 
indicate positive (negative) vortices; (middle): coarse-grained profiles of positive vortex density 
n + (solid red line), negative vortex density n~ (dashed black line) and total vortex density n 
(dot-dashed green line) at t = 0. In the inset, the corresponding coarse-grained profile of the 
polarization magnitude p(y) is reported (solid magenta line); (bottom) coarse-grained profiles of 
superfluid velocity u s (solid red line), normal fluid velocity u n (solid blue line) and counterflow 

velocity u ns = u n — u s (solid green line) at t = 0 
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Figure 4: (Color online), (top): vortex distribution at t — Tf, red empty (black filled) circels 
indicate positive (negative) vortices; (middle): coarse-grained profiles of positive vortex density 
n + (solid red line), negative vortex density n~ (dashed black line) and total vortex density n 
(dot-dashed green line) at t = Tf. In the inset, the corresponding coarse-grained profile of the 
polarization magnitude p(y) is reported (solid magenta line); (bottom) coarse-grained profiles of 
superfluid velocity u s (solid red line), normal fluid velocity u n (solid blue line) and counterflow 
velocity u ns = u n — u s (solid green line) at t = Tf. Red and blue dot-dashed lines indicate the 
initial laminar profiles of the superfluid and the normal fluid, respectively. 
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Figure 5: (Color online), (top): coarse-grained profiles of positive vortex density n + (solid red 
line), negative vortex density n~ (dashed black line) and total vortex density n (dot-dashed green 
line) at t = 6.8 x 10 -3 Ty. In the inset, the corresponding coarse-grained profile of the 
polarization magnitude p(y) is reported (solid magenta line); (bottom) coarse-grained profiles of 
superfluid velocity u s (solid red line), normal fluid velocity u n (solid blue line) and counterflow 
velocity u ns = u n — u s (solid green line) at t = 6.8 x 10~ 3 Tj. Red and blue dot-dashed lines 
indicate the initial laminar profiles of the superfluid and the normal fluid, respectively. 


A. Streamwise flow of expanding vortex-rings 

The two-dimensional vortex-points motion described in Section IIIIB1 can be physically in¬ 
terpreted in three dimensions as an idealized streamwise flow of expanding vortex-rings lying on 
planes perpendicular to and drifting in opposite direction with respect to v(). This vortex-ring 
three-dimensional analogue of the vortex-points motion stems from the vortex points equations 
of motion © and can be clearly discerned if we consider the motion of an anti-vortex pair whose 
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Figure 6 : (Color online). Coarse-grained profile of the longitudinal component of the mutual 
friction force F 1 at different selected times: t — 0 (dot-dashed green line); t — 6.8 x 10 - 3 I> 

(dashed green line); t — Tf (solid green line). 
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initial configuration is symmetrical with respect to the mid plane of the channel and very close to 
the latter, see Fig. 0 

Let r ± (t) = (x±(t), y±(t)) be the trajectories of the positive and negative vortices which consi- 
titute the anti-vortex pair, with initial condition (Vj_, y± ) . The axisymmetric hypothesis imposes 
| y°_ | = \y {) _ = y°, while the proximity to the channel’s mid plane implies t / 0 C 1. According to 
the proposed parallel, the dynamics of this anti-vortex pair corresponds to the three-dimensional 
motion of a very small circular vortex ring centered on the channel’s mid plane and initial radius 
R° = y°. To obtain the typical motion of the pair of anti-vortices (corresponding to the intersec¬ 
tions of the vortex ring with the two-dimensional channel), we average the vortex-points equation 

dY-±- 

of motion © in the streamwise direction and over time deducing the following equation for -— 

at 


dr± 

dt 


= r± = 


( x±(y) 
\y±(y ) 


(1 - a')u s {y) + a'u n (y) 
±a ( u n (y ) -u s (y)) 



(24) 


where the dot operator indicates the time derivative and u ns = u n — u s , to ease notation. In 
this simple axisymmetric anti-vortex pair model, x± = Ur and y ± = R, where Ur and R are 
the averaged vortex-ring streamwise drifting velocity and its expansion rate, respectively. We 
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Figure 7: (Color online). Trajectories r + (t) and r .(t) of an anti-vortex pair, whose initial 
configuration is symmetric with respect to the mid plane of the channel (in dashed blue line). The 
time interval between consecutives positions is constant, indicating an increasing streamwise 
velocity as the vortex-points approach the walls. In an axisymmetric interpretation, the 
three-dimensional analogue of this two-dimensional motion is a streamwise flow of an 

expanding vortex ring. 


therefore have the following relations: 

u R (y) = (1 - a')u a (y) + a'u n (y ) (25) 

R(y) = ±au ns (y ) (26) 

From equation [261 given the plot of u ns (y ) reported in Fig. [4] i.e. u ns (y) < 0 Vy, it clearly 
emerges that the positive (negative) vortex moves towards the y = — 1 (y = 1). Hence, only the 
three-dimensional corresponding vortex rings whose circulation is oriented in the same direction 
of expand, while vortex-rings of opposite circulation always shrink. The trajectory of an 
expanding anti-vortex pair is reported in Fig. [7] 

We would like to stress, however, that our numerical simulations grasp a more general and 
complex dynamics, not enforcing an axisymmetric vortex-points motion, but moving each vor¬ 
tex individually. Therefore, the idealized three-dimensional vortex-ring motion described in the 
present paragraph is a physical interpretation of the average vortex-points motion only. We reckon, 
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nevertheless, that it describes three-dimensionally the most relevant characteristics of the two- 
dimensional flow analyzed in the present model. In Fig. [8] several three-dimensional physical 
interpretations of the vortex-points motion are illustrated, with Fig. [8] (c) describing the vortex 
ring analog of the comprehensive two-dimensional motion described in Section fill Bl 



Figure 8: (Color online). Distinct three-dimensional interpretations of vortex-points motion: (a) 
straight vortices model; (b) idealized axisymmetric vortex ring interpretation described in Section 
IIV At (c) vortex ring analog of the comprehensive two-dimensional motion. 


B. Congruity with vortex-tangle dynamics 

Helium II counterflows are well known to exhibit anisotropic characteristics: the vortex lines 
tend to lie on planes perpendicular to v£. This can be easily deduced, for instance, by the plot 
of the projection of the vortex-line length in the streamwise direction (A x ) in^2 and the plots of 
the anisotropic parameter I' As a consequence, we reckon that our idelized vortex-rings- 
flow model is able to capture the dynamics of the most relevant fraction of the vortex-tangle. In 
addition, it is worth emphasizing that the vortex-lines aligned in the streamwise direction (which 
we neglect in our simplified three-dimensional interpretation) are only affected very slightly by 
the mutual friction interaction which governs the vortex-tangle dynamics. On the other hand, our 
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model is less reliable in the near-wall region where the vortex-tangle assumes a more isotropic 
character. 

Furthermore, from Eqs. (l25l) and (l26l) and the plots of u s , u n and u ns reported in Fig. @]it 
is possible to deduce that in the proposed three-dimensional physical interpretation of our two- 
dimensional model, the vortex-rings drifting velocity in the streamwise direction increases as the 
radius of the vortex-rings grows ( i.e . as the vortex-rings approach the channel walls). This vortex 
dynamics also emerges from past numerical three-dimensional studie d — which describe the 
vortex lines moving towards the solid boundaries with increasing streamwise velocity in opposite 
direction with respect to the normal fluid flow. 

To conclude this section, it is important to underline that the orientation of the expanding 
vortex-rings (circulation in the same direction of the normal fluid mean flow) is responsible for 
the non-uniform profile of u s illustrated in Fig. 01 the superfluid velocity field induced by such 
vortex-rings slows down the superflow in the central region of the channel while the image vor¬ 
tices increase the superfluid velocity near the boundaries. This non-uniform superfluid velocity 
profile is qualitatively recovered in past numerical simulations^. 

Having described what we propose is the three-dimensional physical interpretation of the 
vortex-points motion numerically investigated in our simulations and having discussed its con¬ 
sistency with the vortex-tangle dynamics observed in past three-dimensional numerical studies, 
we reckon that our model, although being two-dimensional, is capable of grasping the most es¬ 
sential and relevant dynamics taking place in helium II T-I channel counterflows. 


V. CONCLUSIONS 

In this work we have performed two-dimensional self-consistent, coupled numerical simula¬ 
tions of helium II channel counterflows with corresponding vortex-line density typical of counter¬ 
flow experiments.-^^ 

The main features of our model are the presence of solid boundaries and the dynamical cou¬ 
pling of vortices and normal fluid. These features make our model more realistic than previous 
investigations, although, due to computational constraints we had to use a two-dimensional ge¬ 
ometry rather than a three-dimensional one. We reckon, however, that our model, despite its 
reduced dimensionality, is capable of grasping, at least qualitatively, the most relevant features of 
the vortex-tangle dynamics occurring in helium II T-I counterflows. For instance, the proposed 
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physical three-dimensional interpretation of the vortex-points motion ( i.e . a streamwise flow 
of expanding vortex rings) is qualitatively in agreement with the three-dimensional vortex-lines 
motion computed under prescribed normal fluid flow4^ In addition, the vortex density profiles 
computed three-dimensionally with imposed Poiseuille normal fluid How— — are consistent with 
the profiles calculated in our two-dimensional simulations. Experimentally, these profiles could be 
estimated by suitable second sound attentuation measurements, employing high harmonics waves. 

In conclusion, the numerical results achieved in our work confirm the already observed velocity 
matching^ and vorticity locking 4142 ’ 44 ’ 63 between the two helium II components. Above all, our 
numerical model predicts the shape of the profile of the normal fluid which has been just observed 
experimentally in channels using laser-induced fluorescence of metastable helium molecules.— 
Furthermore, our results are useful for the interpretation of actual and future experiments, includ¬ 
ing pure supcrllow— and the motion of tracer particles. -^21 
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